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Abstract 

Recently there has been a great deal of interest in studying the properties of dense quark matter, 
with particular reference to diquark condensates and color superconductivity. In the present work 
we report calculations made for the excitations of quark matter for relatively low densities of the 
deconfined phase and in the absence of meson or diquark condensation. Here, we are interested in 
elucidating the role of “Pauli blocking”, as such blocking affects the calculation vector, scalar and 
pseudoscalar qq excitations. As a byproduct of our analysis, we extend our calculations to higher 
densities and explore some consequences of the use of density-dependent coupling parameters for 
the Nambu-Jona-Lasinio model. (Such density-dependent parameters have been used in some of 
our previous work.) For our analysis made at large values of the matter density, we assume that at 
about 13 times nuclear matter density quark matter has only minimal nonperturbative interactions. 
At that high density we compare the result for hadronic current correlation functions calculated 
with density-dependent and density-independent NJL coupling constants. We hnd evidence for the 
use of density-dependent parameters, since the results with the density-independent constants do 
not go over to the perturbative description which we assume to be correct for p ~ Qpc, where pc is 
the matter density for the hnite-density confinement-deconfinement transition. The use of density- 
dependent coupling constants in the study of diquark condensates and color superconductivity has 
not been explored as yet, and is a topic requiring further investigation, particularly given the strong 
interest in the properties of dense quark matter and color superconductivity. 
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I. INTRODUCTION 


At present there is a great deal of interest in exploring the properties of quantum chromo¬ 
dynamics (QCD) at high temperature and baryon density. We may quote Renk, Schneider 
and Weise [1]: “Evidently, investigating the changes of spectral distributions of pseudoscalar 
and vector (as well as axial-vector) excitations of the QCD vacuum with changing tempera¬ 
tures and baryon densities, from moderate to extreme, is a key to understanding QCD ther¬ 
modynamics, its phases and symmetry breaking patterns.” Studies related to that program 
may carried forward using lattice simulations of QCD or by using effective Lagrangians. 

A good deal is known concerning the properties of QCD at hnite temperature. However, 
the properties of QCD at hnite matter density are not well known, since the lattice simula¬ 
tions of QCD at hnite chemical potential are only in a relatively early state of development 
[2-7]. In earlier works we have used a generalized Nambu-Jona-Lasinio (NJL) model to 
study the conhnement-deconhnement transition for light mesons at both hnite temperature 
[8] and hnite density [9]. We have also studied the excitations of the quark-gluon plasma 
by calculating hadronic current correlation functions at hnite temperature [10, 11]. In the 
present work we report upon calculations of such correlators at hnite matter density in the 
deconhned phase. As in the large number of applications of the NJL model and related mod¬ 
els to the study of high-density matter [12-15], we do not consider explicit gluon degrees of 
freedom. Our calculations do not take into account the formation of diquark condensates 
associated with color superconductivity. One reason for not studying diquark condensates 
and related matters at this time is that we have some concerns as to the use of the NJL 
model with constant values of the coupling constants at high density or high temperature. 
It is easier to discuss the matter of temperature-dependent ehective coupling constants [8, 
10] than to discuss density-dependent coupling constants [9], since a good deal is known 
about QCD thermodynamics at hnite temperature. We denote the critical temperature for 
the conhnement-deconhnement transition as Tc and have used Tc = 170 MeV in our work 
[8, 10]. From studies of a pure gluon system in QCD, it is found that at high temperatures 
one approaches a weakly interacting system slowly, with some nonperturbative ehects still 
present at temperatures T = 3Tc or 4Tc [16]. For dehniteness, we assume that we can ne¬ 
glect nonperturbative ehects for T > GT^. For these high temperature we would expect to 
be able to calculate hadronic current correlators using the NJL model and obtain the result 
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expected in lowest-order perturbation theory at energies for which the effective Lagrangian 
may be used. (In that regard, we limit our applications to < 4 GeV^.) However, we 
found in our studies that the calculated hadronic current correlations functions had res¬ 
onant features for T ~ 6Tc, when temperature-independent coupling constants were used 
[11], That was one of a number of reasons that we replaced the temperature-independent 
coupling constants of the NJL model by G{T) = G[1 — O.ITT/TJ. (We remark that other 
dependence on the parameter T/Tc may be assumed, however, we have only explored the 
linear dependence described here.) By analogy, we may introduce G{p) = G[l—{dp/ where 
Pc is the critical density for the confinement-deconfinement transition. In our earlier work 
[9] we took Pc = 2.25p7VM, where pi^M is the density of nuclear matter. Our introduction 
of the parameter (d was not done in a systematic fashion in Ref. [9], but the values used 
were fd ~ 0.18. For the present work, we have used the same value of pc = 2.25pnm and 
have taken (d = 0.17. Therefore, for p > Opc, we have assumed that the system has limited 
nonperturbative features. Since we are not describing diquark condensation in this work, we 
might limit ourselves to consideration of p values larger than, but not too different from pc- 
However, here we will also study the larger values of p < 6pc in order to gain information 
about a possible density-dependence of the NJL coupling constants, in analogy to what was 
done in the case of our finite temperature studies [8, 10, 11]. 

Our calculations are made for zero values of the chemical potential. Thus, we can present 
our results for hadronic current correlators for various values of p/pc- When calculations 
are made in that manner, we can present a particularly transparent discussion of the role of 
“Pauli blocking”. 

Our model for quark matter is that of two ideal Fermi gases of up and down quarks with 
Fermi momentum pp. The density of the quark matter is given by pq = {2Nc/37r‘^)p^, where 
Nc = 3 is the number of colors. That expression may be put in contrast to the expression 
for the density of nuclear matter, p^M = {2/3n‘^)kp, where kp is the Fermi momentum 
of the ideal gas of nucleons. With our form of G(p) we see that, in our model, we have 
assumed that at pq = 5.88pc = 13.2pjYM the system is weakly interacting, since G(p) = 0 
at that density. We may ask if we can apply the NJL model at such high densities. When 
Pq = 13.2p/\iM, we hud pp = 1.64kp = 0.44 GeV. Here we have used kp = 0.268 GeV for 
nuclear matter. Since a typical (sharp) three-momentum cutoff for the standard NJL model 
is 0.631 GeV [11, 17], we see that it is still possible to use the NJL model at the largest 
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density considered here. (It is possible to extend the range of application by using a larger 
momentum cutoff for the NJL model. However, the cutoff is usually hxed by htting the pion 
decay constant, so that the range of variation of the cutoff parameter is limited [17-19].) 

In this work we study the correlators which involve the excitation of states with the 
quantum numbers of the p, tt, /q and rj mesons. The organization of our work is as follows. 
In Sec. II we provide expressions for the hadronic current correlators for the four cases 
considered and stress the role of Pauli blocking in modifying these expressions. In Sec. Ill 
we provide results for the imaginary parts of the correlators and Cp{P'^) for several 

values of p/Pc- In Sec. IV we present the values obtained for the correlators of flavor octet 
currents, CrjSsiP"^) Cf^^ssiP"^)- In Sec. V we compare the results for the correlators at 
p/pc = 5.88 obtained with constant NJL coupling constants with those obtained for the 
density-dependent coupling constants used in this work and in our earlier work [9]. Finally, 
in Sec. VI we present some further discussions and conclusions. In the following discussion 
p will represent the density of quark matter. However, we will sometimes use the notation 
Pq for that quantity. 


II. CALCULATION OF HADRONIC CURRENT CORRELATION FUNCTIONS 


For ease of reference, we present the Lagrangian of our generalized NJL model that 
incorporates a covariant model of conhnement [20-24] 

^ 8 

C = q{i^ - m^)q + '^[{qPqY + 


i=0 

i=0 

+ ^{det[g(l + 75 )g] + det[g(l - 75 )?]} 

pPconf ■ ( 2 - 1 ) 


Here the A*(i = O,--- ,8) are the Gell-Mann matrices, with A'^ = a/2/3 1 , = 

diag(m°,m/m°) is a matrix of current quark masses and Cconf denotes our model of con¬ 
finement. The fourth term on the right is the ’t Hooft interaction. (The effective coupling 
constants in the channels with tt, /o or p quantum numbers are linear combinations of Gs 
and Gd [17].) 
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FIG. 1: (a) The basic vacuum polarization diagram of the NJL model is shown. The lines represent 
a constituent quark and antiquark of mass m{pq). 

(b) The filled region represents a confinement vertex used for calculations made for p < pc [9]. 

In this work we make use of the density-dependent constituent quark masses that were 
calculated in Refs. [9, 25] using a mean-held approximation. The up (or down) quark mass 
mu{p) is rather small for p > Pc and has relatively little effect when calculating the hadronic 
current correlators. The strange quark mass for p ^ Pc is approximately constant with 
ms{p) — 440 MeV [9]. (The difference between the behavior of the up and down quark mass 
and that of the strange quark mass is due to the fact that the quark matter we consider is 
nonstrange.) 

In Fig. la we show the basic vacuum-polarization diagram of the NJL model, where the 
lines denote either constituent quarks or antiquarks. (In Fig. lb we show the introduction 
of a conhnement vertex which is needed in our studies for T < or p < p^. Since we are 
considering the deconhned phase, we can disregard the conhnement vertex for this work.) 
With reference to Fig. 1, we denote the quark momentum as P/2 + k and the antiquark 
momentum as —P/2 + k, and work in the frame where P = 0. 

In forming the imaginary parts of the vacuum polarization functions, the quark and 
antiquark go on mass shell, so that we have = (P°)^/4 — w?{p). It is easy to see that 
there is a minimum value of \k\ = pp, where pp is the Fermi momentum of either the 
up or down ideal quark gases. For smaller values of |/c| the excitation is “blocked” by the 
Pauli Principle. Thus, we see that in calculating the imaginary part of a vacuum polarization 
function, one obtains nonzero values above a minimum value of P°, (P°)min = 4(p|,-|-m^(p)), 
where pp is related to the quark density by the expression pp = /2Nc)pq, given earlier. 

We now consider the imaginary part of the vacuum polarization function corresponding 
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to scalar currents. We find in the case an up (or down) quark is excited 


ImJ^(p2) 


^4vr ^ 


3/2 

- 4(p/ + mi(p))]. 


( 2 . 2 ) 


This expression contains a factor of 2 arising from the flavor trace. Here /c^ = Pg^/4 — 'nr^(p) 
appears in the Gaussian regulator. (We have used Gaussian regulators for most of our 
calculations made using the NJL model. The value of a = 0.605 GeV yields results that 
are similar to those obtained with the sharp three-momentum cutoff parameter A = 0.631 
GeV.) For pseudoscalar mesons, we may use Eq. (2.2) with the phase-factor exponent of 
3/2 replace by 1/2 when the constituent mass is small [10]. The real parts of the vacuum 
polarization functions are obtained by means of a dispersion relation. 

We now consider the calculation of density-dependent hadronic current correlation func¬ 
tions. The general form of the correlator is a transform of a time-ordered product of currents, 


C{P\p) = i 


« T{j{x)j{0) », 


(2.3) 


where the double bracket is a reminder that we are considering the hnite density case. 

For the study of pseudoscalar states, we may consider currents of the form jp^i{x) = 
q{x)i'y 5 X^q{x) where, in the case of the n mesons, i = 1,2, and 3. For the study of 
pseudoscalar-isoscalar mesons, we again introduce jp,i{x) = q{x)Pq{x), but here i = 0 
for the flavor-singlet current and i = 8 for the flavor-octet current. 

In the case of the tt mesons, the correlator may be expressed in terms of the basic vacuum 
polarization function of the NJL model, Jp(P^,p). Thus, 


cAP^p) 


MP^p)j 


1 

G.(p)Jp(P2,p)’ 


(2.4) 


where G^ip) is the coupling constant appropriate for our study of the n mesons. (We have 
found ^7^(0) = 13.49 GeV“^ by htting the pion mass in a calculation made at p = 0.) 

For a study of the correlators related to the p meson, we introduce conserved vector 
currents jfj.,i{x) = g(a;) 7 ^Ajg(a;) with i = 1,2 and 3. In this case we dehne 

J7(P^p) = J,(P^p), (2.5) 


taking into account the fact that the current j^^i{x) is conserved. We may then use the fact 
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that 


MP\p) = 

^ m W + imKp) 

3 An 

- lup\p) 

and write the approximate expression 

‘2 N f p 
ImJ,(p2,p)~-^h - - 






e[P^ - A{pI + ml{p))]{2.T) 


Aml{p) 




e[P^ - A{pI + ml{p))] 


for the vacuum polarization function of the vector-isovector currents. Here k'^ = Pq/A — 
m\{p) appears in the Gaussian regulator. Thus, we dehne 


^ [ l-GyUh(P^) \ • 

where we have suppressed reference to the density dependence of the correlator and the 
vacuum polarization function. We have used Gv{p) = Gy[l — 0.17p/pc] with Gy = 11.46 
GeV“^ for the calculators reported here. 

The calculation of the correlator for scalar-isoscalar states is more complex, since there are 
both flavor-singlet and flavor-octet states to consider. We may dehne polarization functions 
for u, d and s quarks: Ju{P^, p), p) and Js{P^, p)- These functions do not contain the 

factor of 2 arising from the havor trace that was introduced when calculating ImJ^(P^,p), 
ImJ'^(P^,p) and ImJp(P^,p) earlier in this section. 

In terms of these polarization functions we may then dehne 


= -^[Juip^^p) + Jd{p^,p) + Js{p^,p)\ , 


( 2 . 11 ) 


Jo%{p^,p) = ^[Ju{p^,p) + Jd{p‘^,p) - 2J,(p^p)], 


( 2 . 12 ) 


J8s{p‘^,p) = ^[Ju{p‘^,p) + Jd{p‘^,p) +AJs{p^,p)] . 


We also introduce the matrices 


2 ^ ’Joo{p‘^,p) Jo8{p‘^,p) 


JiP :P) = 


J80{p^,p) ^88(p^p) 


(2.13) 


(2.14) 
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(2.15) 


G(P) 


Goo(p) Go&{.p) 
Gsq{p) Gssip) 


G{p\p) = 


Goo{p^,p) Go8{p^,p) 
Gso{p^,p) Gs8{p^,p) 


and write the matrix relation 


Gip\p) = Jip\p)[l-G{p)J{p‘^,p)] ^ 


(2.16) 


(2.17) 


III. RESULTS OF NUMERICAL CALCULATIONS OF PSEUDOSCALAR AND 
VECTOR HADRONIC CURRENT CORRELATION FUNCTIONS, AND 

Cp{P‘^) 


In Fig. 2 we present values calculated for ImU,r(P^) for p/pc = 1.2, 2.0, 3.0, 4.0 and 5.88. 
The threshold for each curve is given by (P°)mm = ‘^{Pf + ^uiPq))- If Pc = ‘^■‘^^Pnm, we 
have 


Pf 


/2.25 p 

V Nc Pc 


1/3 

Pf , 


(3.1) 


where Pf is the Fermi momentum of nuclear matter, Pp = 0.286 GeV. Thus, when p/Pc = 
5.88, we hnd pp = 0.439 GeV, which is less than the standard three-dimensional (sharp) 
cutoff of A = 0.631 GeV that is often used for the NJL model. It is seen from the hgure 
that there are significant nonperturbative effects, except at p/pc = 5.88, where ImG,r(P^) = 
ImJ.,r(-P^). The last result follows, since G^^i^p) = 0 for p/pc = 5.88. 

In Fig. 3 we show the results of a similar calculation of ImGp(P^). Here the thresholds 
for the various curves are the same as those seen in Fig. 2. 


IV. RESULTS OF NUMERICAL CALCULATION OF SCALAR AND PSEUDO¬ 
SCALAR CORRELATORS G/o88(P^) AND C^88{P^) 

In the case of the correlator of pseudoscalar-isoscalar currents, the components 
Gr/osiJ^^) = Gj^8o(-P^) and Gr,88{P‘^) are all important. We choose to show ImG,,88(-P^) in Fig- 
4 for the values of pjpc = 1.5, 2.0, 3.0, 4.0 and 5.88. The resonant structures have rather 
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FIG. 2: The figure presents values of ImC' 7 r(P^) for various values of p/Pc- Here, p/pc = 1.2 [solid 
line], 2.0 [dashed line], 3.0 [dotted line], 4.0 [dashed-dotted line] and 5.88 [dashed-(double)dotted 
line]. We have used G-,^ = 13.51 GeV“^. 



FIG. 3: The figure shown the values of ImC'p(P^). [See the caption of Fig. 2.] Here we have used 
Gv = 11.46 GeV-2. 








FIG. 4; The figure shows the values calculated for the imaginary part of the correlator of pseu¬ 
doscalar, flavor-octet currents ImC'r;88(P^). Here, pjpc = 1-5 [solid line], 2.0 [dashed line], 3.0 
[dotted line], 4.0 [dashed-dotted line], 5.0 [dashed-(double)dotted line], and 5.88 [short dashed 
line]. 

small widths when compared to the features seen in Fig. 2 and 3. It is also worth noting 
that each curve has two thresholds, one corresponding to m„(p) and the other corresponding 
to m,(p). Thus, we have (P°)Ln = ^\Pf + (^°)Ln = ^[Pf + mlip)]- 

Similar remarks pertain for the scalar current correlators, C/oOo(-P^), C'/o 08 (-P^) = 
C/q8o(-P^) and The values of ImC/(,88(-P^) are shown in Fig. 5 for various values 

of p/Pc- Note that the resonances seen in that hgure are quite narrow. That is probably 
due to the different phase space behavior for the /o (p-wave) and p (s-wave). 


V. RESULTS OF NUMERICAL CALCULATIONS WITH DENSITY- 
DEPENDENT AND DENSITY-INDEPENDENT NJL COUPLING CONSTANTS 

In this section we are interested in presenting some evidence for the density-dependent 
coupling constants used in our work [9]. (As noted earlier, the argument is more easily made 
when we introduce temperature-dependent coupling constants, since much more is known 
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FIG. 5; The values of the imaginary part of the correlator of scalar, flavor-octet hadronic currents, 
ImC'/Q88(.P^) is shown. [See caption of Fig. 4.] 

concerning QCD thermodynamics at finite temperature than at finite density.) In the case 
of finite density, we have assumed that the system is weakly interacting at p = 5.88pc with 
Pc = 2.25pArM- (Other values for pc could be used, but here we will continue to explore the 
consequences of the choice made in our earlier work [9].) In Fig. 6 we compare the results 
of our model for ImO^(P^) at p/pc = 5.88 with the results obtained when we use a constant 
value for As seen in the hgure, there is about a factor of 3 difference in the result of the 
two calculations. The difference in the two calculated results for ImGp(P^) seen in Fig. 7 is 
not quite as marked as that seen in Fig. 6, but is still signihcant. 

In Fig. 8 we compare the results for ImG,j88(-P^) at p/pc = 5.88 for the two methods of 
calculation. A rather dramatic difference is seen in Fig. 8, where we see that the calculation 
with constant values of G^, G^g and G^ leads to a resonance at P^ = 1.53 GeV^. In Fig. 
9 we show similar results for ImG/(,88(-P^) with a resonance at P^ = 1.58 GeV^ in the case 
that constant values of G^, G^ and G|^ are used. 
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FIG. 6: The solid line shows the values of ImC' 7 r(P^) for the case of density-independent NJL cou¬ 
pling constants. The dashed-(double)dotted line represents ImC' 7 r(T’^) calculated with the density- 
dependent coupling constants Gt^{p) = G^[l — O.lTp/pc] for pjpc = 5.88. 

VI. DISCUSSION 

In this work we have assumed that deconfinement takes place at pc = 2.25pnm and that 
the quark gluon plasma is weakly interacting at p = 5.88pc = 13.2pArM- This particular 
model was explored in Ref. [9] where we calculated the masses of several mesons and their 
radial excitations for various matter densities, p < pc, making use of our generalized NJL 
model with conhnement. In the present work we have put forth some evidence that the 
NJL coupling constants should be density dependent to obtain a consistent formalism. We 
may consider what results would be obtained if we used a larger value of the density for 
the conhnement-deconhnement transition. We have investigated the values of pc = S.Opatm, 
Pc = ‘i.QpNM and Pc = S.OpATM and in each case have assumed that the quark-gluon plasma 
is weakly interacting for p = 6.0pc. (We should note that, if pc is made large, the value 
calculated for pp will become larger than the value of the three-momentum NJL cutoff, 
A = 0.631 GeV, that is often used. However, our computer code still provides results under 
that circumstance, since we use a Gaussian regulator.) It is found in all the calculations 


12 






P" (GeV") 


FIG. 7; The solid line shows the values of ImC'p(P^) for the case of density-independent NJL cou¬ 
pling constants. The dashed-(double)dotted line represents ImC'p(P^) calculated with the density- 
dependent coupling constants G'y(p) = Gy)! — O.lT/o/pc] for p/Pc = 5.88. 

made at the values of p for which we have assumed the system to be weakly interacting, 
that the values of the correlators calculated with constant values of the coupling constants 
differ signihcantly from the results obtained with the density-dependent coupling constants. 
We suggest that this matter should be resolved before we undertake studies of diquark 
condensation at high densities and related matters. 

One characteristic of our calculation of hadronic current correlators at hnite temperature 
and hnite density in the deconhned phase is the appearance of complex resonance structure 
in some cases. That feature is in general agreement with the observation made in Ref. 
[26] “...that the correlators possess a nontrivial structure in the deconhned phase.” Meson 
correlators were studied in hnite temperature lattice QCD in Ref. [27], in which the authors 
state, “...Below Tc we observe little change in the meson properties as compared with T = 
0. Above Tc we observe new features: chiral symmetry restoration and signals of plasma 
formation, but also an indication of persisting “mesonic” (metastable) states and diherent 
temporal and spacial masses in the mesonic channels. This snggests a complex pictnre of 
QGP in the region (l-1.5)Tc”. 
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P" (GeV') 


FIG. 8; The figure shows ImC'^88(-P^) calculated for p/pc = 5.88. The dashed-(double)dotted line 
is the result of using our density-dependent coupling parameters, while the solid line is the result 
when density-independent coupling parameters are used. 



P" (GeV") 


FIG. 9: The figure shows ImC'/p88(F’^) 
calculated for pjpc = 5.88. [See the caption of Fig. 8.] 
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